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-»-» ■ ABSTRACT 

This paper exploits the gravitational magnification of SNe la to measure properties 
of dark matter haloes. Gravitationally magnified and de-magnified SNe la should be 
brighter and fainter than average, respectively. The magnification of individual SNe la 
can be computed using observed properties of foreground galaxies and dark matter 
halo models. We model the dark matter haloes of the galaxies as truncated singular 
isothermal spheres with velocity dispersion and truncation radius obeying luminosity 
dependent scaling laws. 

A homogeneously selected sample of 175 SNe la from the first 3-years of the Su- 
pernova Legacy Survey (SNLS) in the redshift range 0.2 < z < 1 is used to constrain 
models of the dark matter haloes associated with foreground galaxies. The best-fitting 
velocity dispersion scaling law agrees well with galaxy-galaxy lensing measurements. 
We further find that the normalisation of the velocity dispersion of passive and star 
forming galaxies are consistent with empirical Faber-Jackson and Tully-Fisher rela- 
tions, respectively. If we make no assumption on the normalisation of these relations, 
k>( \ we find that the data prefer gravitational lensing at the 92 per cent confidence level. 

• Using recent models of dust extinction we deduce that the impact of this effect on our 

results is very small. 

We also investigate the brightness scatter of SNe la due to gravitational lensing, 
which has implications for SN la cosmology. The gravitational lensing scatter is ap- 
proximately proportional to the SN la redshift. We find the constant of proportionality 
to be B ~ 0.055l°;°4? ma S (B ^ 0-12 mag at the 95 per cent confidence level). If this 
model is correct, the contribution from lensing to the intrinsic brightness scatter of 
SNe la is small for the SNLS sample. According to the best-fitting dark matter model, 
gravitational lensing should, however, contribute significantly to the brightness scatter 
at z> 1.6. 

Key words: gravitational lensing - supernovae: general - dark matter - galaxies: 
haloes 
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1 INTRODUCTION 

In the context of supernova cosmology, gravitational lensing 
is usually regarded as a source of uncertainty because 
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it adds extra scatte r to the brightness of high redshift 
Type la s upern o vae dKantowski. V aughan & Branch 19951: 
Friemanl Il997l; IWambsganss et all Il997l; Itiolz fc Waldl 



19981: iBergstrom et all l2000l : IhoIz fc Lindeil 120051 : 
Gunnarsson et al. I l2006t ). Due to flux conservation, the 
effects of gravitational lensing magnification and de- 
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magnification average out and is therefore expected to iead 
to n egligible bias in cosmological parameter estimation (see, 
e.g., ISarkar et al1l2008l : I Jonsson et ai1l2008h . 

Gravitationally lensed Type la supernovae (SNe la) 
are, however, interesting in their own right, because they 
can be used to measure magnification. Magnified and de- 
magnified standard candles, like SNe la calibrated using 
light curve shape and colour corrections, should appear to 
be brighter and fainter than average, respe ctively. Cosmic 
magni fication has already been detected bv lScranton et all 
(2005) using the correlation between quasars and galax- 
ies observed within the Sloan Digital Sky Survey. The 
magnification (or de- magnification) of SNe la can be used 
to probe cosmology ( MetcalJ Il999l; iDodelson fc Vallinottol 



120061 ; IZentner fc Bhattach arva 20 09|) and the nature of the 
lensing matter (lRauchlll99l|; iMetcalJ Il998l; fMetcalf fc Silk] 
1999 1; ISeliak fc Hold Il999l; llVIortsell. Goobar fc Bergstronl 
200ll ; lMintv et al.ll2002l ; iMetcalf fc SilkH2007l ; I Jonsson et all 
2009h . 

In addition to this paper there is another l ensing study 
maki ng use of Supernova Legacy Survey (SNLS, Astier et al.l 
I2006T ) data. The study of lKronborg et all l|2010l ). which fo- 
cus on the detection of a gravitational lensing signal using 
assumptions about the dark matter halo models, report the 
detection of a signal at the 99 per cent confidence level. This 
result confirms an earlier tentative detection with GOODS 
data l| Jonsson et al.ll2007f ). In the work presented here the 
assumptions of the underlying model are tested. 

We use presumably lensed high-redshift (0.2 < z < 1) 
SNe la from the SNLS, to measure properties of dark mat- 
ter haloes of galaxies in the deep Canada-France-Hawaii 
Telescope Legacy Survey (CFH TLS) fie l ds. T he method 
we use was f i rst su ggested by iMetcalll (| 19981 ). Recently 
IJonsson et alj (|2009T ) applied this method to SNe la and 
galaxies observed within the Great Observatories Origins 
Deep Survey (GOODS). Since the SNLS SNe la are not as 
distant (z < 1) as the GOODS SNe la, the effect of gravita- 
tional lensing is expected to be smaller than for the GOODS 
sample. However the SNLS SNe la are far more numerous 
and selected in a more homogeneous way. 

The data used to put constraints on dark matter halo 
properties are described in section[2] In section[3]the method 
is described. The results are presented in section [3] We in- 
vestigate the effect of dust extinction in section [S] Section [5] 
is devoted to gravitational lensing brightness scatter. The 
results are summarised and discussed in section [7] 

Throughout the paper we assume a cosmological model 
characterised by Qm = 0.27, Qa = 0.73, and h = 0.7. 



measure their redshifts ll Howell et al.l 120051 ; iBronder et alj 
120081 ; iBalland et alj I2009T 1. To measure properties of dark 
matter galaxy haloes using lensed SNe la, observations of 
both supernovae and galaxies are required. 



2.1 Galaxies 

The number of galaxies in each field is of the order 10 5 
and we therefore have to rely upon photometric redshifts. 
The techniques us ed to obtain galaxy properties are de- 
scribed in detail in lSullivan et al. (|2006l). Galaxy photome - 
try is per formed using S E xtract or ( Bertin fc Arnou"tslll996l ) 
following ISullivan et all |2006l ). We estimate the galaxy 
redshift, z ga i, and its physical properties using the pho- 
tometric redshift code Z-PEG, based on the PEGASE.2 
galax y spectral evolution code (jFioc fc Rocca-Volmerangel 
119971 ; iLe Borgne fc Rocca-Volmerangel |2002| ). The best- 
fitting spectral energy distribution (SED) is determined us- 
ing a x 2 minimisation between the observed fluxes, the cor- 
responding flux errors, and the synthetic photometry gen- 
erated by integrating the template SEDs through the SNLS 
effective filter responses. 

We then use this best-fitting SED for the estimation 
of the various properties which can be used to characterise 
the galaxy. We first require an appropriate luminosity. The 
galaxies could, in principle, be characterised by the lumi- 
nosity in any rest frame wavelength band, but to facilitate 
comparison with results in the literature we use Mb, the 
absolute magnitude in the B-band. Given the photometric 
redshift of a galaxy and the best-fitting SED, we compute 
absolute magnitudes by integrating the de-redshifted SED 
through the B-band filter. 

We also require a galaxy type r as we expect the rela- 
tions between dark matter halo es and galaxies to be different 
for different galaxy types. As in lSullivan et al.l (|2006l ). we use 
the specific star- formation rate (sSFR), the star-formation 
rate (in Mq yr _1 ) per unit stellar mass. If log(sSFR) < —12, 
the galaxy is classified as passive, otherwise it is classified 
as star forming. 

According to lllbert et al.l ()2006h , reliable photometric 
redshifts can only be obtained for galaxies with m, < 24 
mag from the CFHTLS deep data. We therefore include only 
galaxies brighter than this magnitude limit. The fraction of 
passive galaxies ranges from 5 per cent in Dl and D3 to 
8 per cent in D4. We have checked that including galaxies 
fainter than 7?v = 24 mag, which have less secure redshifts, 
have negligible impact on the results. 



2 DATA 

The SNLS is made up of photometric and spectroscopic ob- 
servations. The photometry is obtained as part of the deep 
component of CFHTLS with th e one square degree imager 
MegaCam (jBoulade et al.ll2003T ). The deep part of CFHTLS 
comprises four fields (Dl, D2, D3, and D4), each ~ 1 deg 2 
in size, imaged in u* , g' , r' , i' , and z' filters approximately 
every 4-5 days during dark and grey time, sui table for de- 
tecting supernovae an d building light curves (|Astier et al.l 
120061 ; iGuv et al] 120091 ) . The spectroscopic observations are 
used to determine the nature of supernova candidates and 



2.2 Supernovae 

We assume that the brightness of a SN la traces its magnifi- 
cation. In our analysis we use the Hubble diagram residual, 

AmsN = (J-b — h(Qm, ^a; zsn), (1) 

and its uncertainty, o"Am SN , to characterise the brightness of 
a SN la. In the equation above, /ib is a distance indicator 
and h(Qm, Qa; zsn) is the predicted value for a SN la at 
redshift «sn assuming a cosmological model with Qm = 0.27 
and Qa = 0.73. The Hubble diagram residual measures the 
brightness of an individual SN relative to the average of the 
population - brighter SNe have more negative residuals. 
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Table 1. Number of SNe la in the CFHTLS deep fields. 



Field 


All SNe la 


Selected SNe la 


Dl 


55 


42 


D2 


19 


33 


D3 


75 


55 


D4 


64 


45 


Total 


243 


175 



Our set of resi duals was computed using the SiFTO 
i|Conlev et al l |200S|) li grit curve fitting package. We have 
also performed the analysis wit h residuals computed us- 
ing the SALT2 (|Guv et all 120071 ') package. The results ob- 
tained using SiFTO and SALT2 residuals are very similar 
and we therefore only present results obtained using SiFTO 
residuals here. For this set of residuals the intrinsic scatter, 
i.e. the component of the scatter which cannot be attributed 
to known observ ational errors, is approximately 0.09 mag 
l|Guv et al.ll2009T l. 

In order to compute the magnification due to the galax- 
ies in the foreground we also need, in addition to the red- 
shift, zsn, the position on the sky, 9sn, for each SN la. When 
available we use the spectroscopic redshift of the host galaxy, 
otherwise the spectroscopic redshift of the SN la itself. Posi- 
tions of the SNe la, as well as of the galaxies, were obtained 
from the CFHTLS i'-band photometry. 

We use 243 SNe la from the 3-year SNLS data set 
|Cuv et al.ll2009h . Since some parts of the deep CFHTLS 
fields are covered by bright stars, they have to be masked 
l|Sullivan et al.ll2006j ). SNe la located too close to a masked 
region or the boundary of the field are removed from the 
sample because of the lack of observations of foreground 
galaxies. To characterise th e size of a masked region we use 
an effective radius, 9 e g = A/n, where A is the area of the 
masked region. The effective radius ranges from a few to a 
few hundred arc seconds. Roughly 20 per cent of the deep 
fields are covered by masked regions, but most of this area 
is made up of a few regions with 9 c a > 60". 

When computing the magnification, all galaxies within 
a radius 9 C a re included. We use 9 C = 60", which according t o 
simulations ijGunnarsson et al ] |2006l : iKronborg et"ai1l2010h , 
is sufficient. For a SN la to be included in the analysis we 
therefore use the selection criterion that the distance to the 
boundary of the field or the distance to the boundary of 
the nearest masked region must be larger than 9 C . If all 
masked regions were taken into account, most SNe la would 
lie too close to a masked region and for that reason fail the 
selection criterion. We therefore ignore all masked regions 
with 6> e ff < 20". Since the location of the masked regions 
should be uncorrelated with the location of the foreground 
galaxies, ignoring small masked regions should not bias the 
results, only increase the uncertainty. 

Only 175 out of the 243 SNe la fulfill the selection cri- 
terion. In Table [1] the total number of SNe la together with 
the number of selected ones are listed for each field. 



puted for a range of halo models. Then the best-fitting halo 
model to the Hubble diagram residuals is obtained by \ 2 
minimisation. In the following the method is described in 
more detail. 



3.1 Computing the magnification 

In order to constrain properties of dark matter haloes we 
exploit the correlation between the brightness of SNe la 
and their magnification. For each SN la the gravita- 
tional magnification due to the dark matter haloes as- 
sociated with the observed galaxies in the foreground, 
Amiens(AsN, Ai oa ; Ahaio), is computed assuming a halo model 
with parameters described by Ahaio- To describe the redshift, 
zsn, and location, #sn, of the SN la, which are used in the 
calculation, we use the shorthand notation 

Asn = {zsn, 0sn}- (2) 

For the calculation of the magnification we also need proper- 
ties of the A^gai galaxies along the SN la line of sight, which 
we denote by 

A, os = {4„ 0U, M%,t\ zy\B^\M^\r N ^}, (3) 

where Zg al , # ga i, M B , and r % are redshift, position on the sky, 
absolute S-band magnitude, and type of the ith galaxy. 

To compute the magnification of a SN la we 
use the weak lensing approxi mation (see, for example, 
ISchneider. Ehlers fc Falcall992l 'l. We have checked the va- 
lidity of this approximation against a more time consuming 
ray tracing algorithm for a few cases. For the sample of 
SNe la considered here, which have relatively low redshifts, 
the difference is found to be < 5 per cent. 

According to the weak lensing approximation the mag- 
nification, in terms of magnitudes, is related to the conver- 
gence, k, (which is a dimensionless surface density) via 

Amiens ^ -2.17k, (4) 

where 

K = Klos ^b- (5) 

The first term in equation ([5| is a sum over the contribution 
to the convergence from each galaxy, K ga i, along the line of 
sight, 

K los = ^ K gal- (6) 
i = l 

The second term, representing the compensating effect of 
the background density, acts as a normalisation allowing the 
magnification relative to a homogeneous universe to be com- 
puted. Due to flux conservation (k) = 0. To ensure that this 
condition is fulfilled we use Kh = (ftios). The average value 
of (tios is computed for a large number of randomly selected 
lines of sight. One important element of the model remains 
to specify, namely K ga i. 



3 METHOD 

The method (|MetcalJ [l998) we use can be summarised in 
two steps. First the magnifications of the SNe la are com- 



3.2 Halo models 

The contribution to the dimensionless convergence from 
each galaxy is given by 



4 J. Jonsson et al. 



K gal(0 = 



(7) 



where the surface density, £(£), is obtained by projecting 
the matter distribution onto a lens plane, 



p{i,y)dy, 



(8) 



where £ is a vector in the plane and y is a coordinate along 
the line of sight. This is the step in the calculations where the 
density profile, p(r), of the dark matter halo enters. Equa- 
tion (0 depends also on the critical surface density, 

^ - AitG D d D ds ' W 

which, in turn, depends on the angular diameter distances 
between the observer and the source, D s , the observer and 
the lens (deflector), Dd, and the lens and the source, Dds- 
The distances are computed from zsn and z m \ assuming a 
cosmological model (fijyr = 0.27, J1a = 0.73, and ft = 0.7). 

We assume that the dark matter haloes can be described 
by a singular isothermal sphere (SIS). The density profile of 
a SIS, 



p(r) 



(10) 



2txG r 2 ' 

has only one free parameter - the velocity dispersion, a. The 
formula for the convergence of a SIS can be considerably 
simplified, if we introduce a suitable length scale £o in the 
lens plane, 



£ = 47TO" 



D s 



(11) 



Since the SIS profile is divergent, the profile is truncated at 
a radius rt- The convergence of a truncated SIS is given by 



Kgal(lE) 







• arctan ^Jx 2 /x 2 



if x ^ x t 
if x > x t , 



(12) 



where x = £/£ and xt = r t /£o- 

To investigate the relationship between galaxy luminos- 
ity and velocity dispersion, we use the scaling law 



(13) 



where L* is a fiducial luminosity, which we take to be 
L* = 10 10 ft" 2 Lq in the S-band. In terms of absolute B- 
band magnitudes, which we will work with, the scaling re- 
lation becomes 



(J — (7*10 



- n (Mg-M* B )/2.5 



(14) 



where M* B = -19.52 + 51og 10 h. 

For the truncation radius we will consider a scaling law 
of the form 

L \ vl 



n = r* ( — ) = r 



(15) 



Since we only include galaxies located a distance 9 C from the 
position of the SN la, the truncation radius has an effect only 
if r t /Dd < 6c- For C — 60" and z ga i < 1 this corresponds to 
r t < 300/i" 1 kpc. 



3.3 Finding the best-fitting model 

In summary, the parameters of our halo model are 



Ahaio = {cr*,r)> r *,l}, 



(16) 



where {a,,ri} and {r*,7} describe the velocity dispersion 
and truncation radius scaling law, respectively. To find 
the best-fitting halo model we minimise the following \ 2 ~ 
function 



[Am SN - Amiens (Ag 



2 \ ^ I DIN "DIN 5 " "lOSl " "nciiu/ I 

X a 2 ) 2 ' (17) 

where the superscript i refers to the ith SN la. 

In equation (|17p the intrinsic dispersion, <7i nt , is added 
in quadrature to the uncertainty of the Hubble diagram 
residual. The uncertainty in Ami ons (AsN, Ai os ; Ahaio), which 
for a specific halo model Ahaio only depends on the uncer- 
tainties in Asn and Ai os , is negligible compared to (j- 1tA and 
(7Am Slg and is for that reason neglected in equation (|17[) . 

3.4 Cosmic variance 

The four deep CFHTLS fields are located widely separated 
in the sky. Due to cosmic variance we expect the density 
of galaxies to vary between the fields. This variation affects 
the normalisation of the convergence, Kb- Figure [T] shows Kb 
as a function of source redshift for the four fields computed 
for two different SIS models. For both models rj = 1/3 and 
7 = 0. The lower set of curves correspond to a* = 130 
km s -1 and r* = AOhT 1 kpc, while the upper set of curves 
corresponds to er„ = 130 km s -1 and r* = 200ft -1 kpc. 
For both models the differences between the fields increase 
with source redshift. In order to reduce the effect of cosmic 
variance, we use the average value of Mifj . £LS indicated by the 
circles 



4 RESULTS 

The halo model, Ahaio, has four parameters. Instead of con- 
straining all parameters simultaneously, which would be 
quite time consuming, we focus on different aspects of the 
model in the following. 



4.1 Truncation radius 

Our model of the luminosity dependent truncation radius, 
equation (|15p . has two parameters: r* and 7. In order to 
reduce the size of the parameter space to a manageable size, 
the investigation is limited to a single value of 77 in equa- 
tion (|13[1 . We expect this parameter to be in the range 1/4 < 
rj < 1/3, since r\ ~ 1/4 and r\ ~ 1/ 3 are the exponents be - 
longing to the Faber-J ackson (F-J, Faber fc Jackson! 0*9761 ) 
and Tully-Fisher (T-F. lTullv fc Fisherlll977h relations valid 
for early and late type galaxies, respectively. Since most of 
the galaxies in our sample are star forming and therefore 
should be of late type, we assume rj = 1/3. We have checked 
that the dependence of rj is rather weak, justifying our re- 
stricted investigation. 

When minimising equation (|17|l with respect to the 
three parameters we find two minima with similar x 2 . For 
the first minimum [a = 114 km s -1 , r* = 170ft -1 kpc, 
7 = -2.8) the value of the \ 2 is 164.50. The \ 2 for the sec- 
ond minimum (a — 124 km s , r* = 69ft _ kpc, 7 = —0.9) 
is slightly higher, 164.66. The best-fitting velocity dispersion 
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Figure 1. Convergence normalisation, function of source 

redshift. Solid, dashed, dash-dotted, and dotted curves corre- 
spond to Dl, D2, D3, and D4, respectively. The two sets of 
curves are computed for two different SIS models. For both models 
r\ = 1/3 and 7 = 0, but the values of a, (measured in km s -1 ) and 
r» (measured in h^ 1 kpc) differ. The lower and upper set of curves 
correspond to {cr,,r»} = {130,200} and {a,,r,} = {130,40}, re- 
spectively. Open circles indicate the average of kj, computed for 
the four fields. 



normalisations, cr», are similar for the two minima. Nega- 
tive values of 7 correspond to rather peculiar truncation 
radius scaling laws. For these scaling laws more luminous, 
and presumably more massive, galaxies have smaller trun- 
cation radii than less luminous galaxies. Unfortunately, the 
data do not allow us to put any limits on r* and 7. 

Fortunately, limits can be obtained for the normalisa- 
tion of the velocity dispersion. Figure [5] shows the probabil- 
ity density function of <r» after marginalisation over r* and 
7. We find <r» = 83lgg km s -1 (68.3 per cent confidence 
level) and er* < 147 km s -1 at 95 per cent confidence level. 

Let us now compare the best-fitting values of the trun- 
cation radius with the virial radius. The virial radius is de- 
fined as the radius inside which the mean density of the 
halo is A v i r (z) times the mean matter density, pm(z) = 



^mPc(1 + 2 ) 3 , where p° is the present critical density. To 
compute the viria l over d ensity we use the approximation of 
iBrvan fc Norman! j 19981 ). 



A vir (z) ~ (18tt 2 + 82x - 39x 2 )p c (z)/p m (z), 



(18) 



where x = p m (z)/p c (z) 
is 



- 1. For a SIS profile the virial radius 



= 20cr [A vil {z)tl M {l + zf] 1/2 ft -1 kpc 



(19) 



From equations (0 and ([9]) it is clear that there is an op- 
timum lens redshift, z ga i, for every source redshift, zsn, de- 
pending on the distance combination D^D^ S / D s , where the 
convergence is maximal. Our distribution of sources peaks at 
zsn — 0.7, which corresponds to an optimum lens redshift of 



b 




250 



Figure 2. Probability density function of <r* after marginalisation 
over r* and 7. The exponent of the luminosity dependent velocity 
dispersion scaling law is rj = 1/3. 



Zgai ~ 0.3. For a dark matter halo at this redshift, the virial 
radius is, according to equation (|19[) , »Vir — 1.6<rft -1 kpc. 

The virial radius of an L* galaxy, for which we find <r» ~ 
120 km s , is hence r* ir ~ 190ft -1 kpc. When exploring 
the truncation radius scaling law we find r* = 170ft -1 kpc, 
which is similar to the virial radius, but the exponent is 
peculiar, 7 = —2.8. As we will see in the next section, we 
find even smaller values of r* for 7 = 0. 



4.2 Velocity dispersion scaling law 

Let us now investigate the velocity dispersion scaling law, 
which is described by the parameters a* and rj. Here we 
consider a constant truncation radius, r t = r* , i.e. we assume 
7 = 0. The best-fitting values of the three parameters to 
the data are a, = 122 km s -1 , 77 = 0.45, and r, = 40ft -1 
kpc. This value of r» is only one fifth of the virial radius 
corresponding to <t*. For this minimum the \ 2 is 164.30. In 
order to obtain constraints on the velocity dispersion scaling 
law (cr* and 77), we marginalise over r*. Reasonable fits to 
the data can not be obtained for all values of r» and we 
therefore apply a flat prior, 



7r(r«) 



if 20ft -1 < r* < 100ft -1 
otherwise, 



(20) 



when marginalising over this parameter. Figure [3] shows the 
results for the velocity dispersion scaling law. The circle indi- 
cates the best-fitting values (c* = 112 km s -1 and r\ — 0.45). 
In the figure are also shown confidence level contours at the 
68.3, 95, and 99 per cent confidence level. 

Figure [4] shows the probability density function of rj af- 
ter marginalisation over cr„ and r*. We find 77 = 0.45:^2? 
(68.3 per cent confidence level) and —0.34 < 77 < 1.17 at 
the 95 per cent confidence level. For the velocity disper- 
sion normalisation we can only obtain upper limits. After 
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Figure 3. Best-fitting values (circle) and 68.3, 95, and 99 per cent 
confidence level contours for the luminosity dependent velocity 
dispersion scaling law after marginalisation over r* . 
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Figure 4. Probability density function of rj after marginalisation 
over (j * and r*. 



marginalisation over r\ and r* we find <r» < 156 at the 95 
per cent confidence level. 

We have redone the fit of the velocity dispersion scaling 
law assuming the truncation radius to be equal to the virial 
radius, given by equation (|19|) . The best-fitting model (x 2 = 
168.13) is in this case cr* = 56 km s _1 and r] — 0.25. After 
marginalisation we find a, < 91 km s _1 at the 95 per cent 
confidence level and r) = Q.2bt° Q 2 2 t (-0.27 < r\ < 1.24 at the 
95 per cent confidence level). 



4.3 Disentangling the velocity dispersion of 
different galaxy types 

As already noted we expect different scaling laws for differ- 
ent galaxy types. The galaxies in our sample are classified 
as either passive or star forming depending on their sSFR 
(see section fO} . We will associate passive and star form- 
ing galaxies with early and late type galaxies, respectively. 
For passive and star forming galaxies we assume the velocity 
dispersion to scale with luminosity as 



and 



sf sf ( L 

a = fJ * \T. 



1/4 



1/3 



(21) 



(22) 



respectively. The exponents of these scaling laws correspond 
to the F-J and T-F relations. Once again we restrict our 
study to models with 7 = 0. 

For the best-fitting halo model, al = 196 km s _1 , 
of = 109 kms" 1 , and r* = 45ft" 1 kpc, the \ 2 is 164.59. Fig- 
ure [5] shows a magnification- residual diagram computed for 
this model. We expect the points to scatter around the solid 
line with slope unity. The scatter around this line is 0.147 
mag and the linear correlation coefficient IS fcorr = 0.22. For 
the Hubble diagram residuals the scatter is 0.150 mag. Cor- 
recting this sample of SNe la for lensing l|Gunnarsson et al.l 



120061 ; I Jonsson. Mortsell fc Sollermanl I2009T ) using the best- 
fitting model would hence lead to a very small reduction of 
the scatter. 

We marginalise over r* using the prior n(r*) to find con- 
straints on al and of. Figure [6] shows the result of our effort 
to disentangle the velocity dispersion for different galaxy 
types. The best-fitting values (erf = 181 and al = 87 km 
s _1 ) after marginalisation over r* are indicated by the circle. 
The point a% — of = 0, corresponding to no lensing, is ex- 
cluded at the 91.6 per cent confidence level. This detection i s 
in good agreement with the results of Ijonsson" et all i|2008r i. 
where the probability of detecting lensing at this confidence 
level for a sample of 175 SNLS SNe la was predicted to be 
80-90 per cent. 

From Fig. [6] it is clear that larger values of the velocity 
dispersion are admitted for passive than for star forming 
galaxies. 

If we marginalise over r* and of {al), we can obtain 



We find al 



91 H 



km 



f st _ ni+36\ 
86 \°* — y± -63/' 

3.3 per cent confidence level) and al < 260 (<7» f < 155) 
1 at the 95 per cent confidence level. 



limits on ctJ (of 



s 

km : 



When assuming r\ — r v i r we find the best-fitting ve- 
locity dispersion normalisation for passive and star forming 
galaxies (x 2 = 167.49) to be al = 113 and af — 33 km 
s~ , respectively. After marginalisation we find al = 1091*0 
km s _1 (al < 172 km s _1 at the 95 per cent confidence 
level) for passive galaxies. For star forming galaxies we ob- 
tain an upper limit at the 95 per cent confidence level after 
marginalisation, cr» f < 101 km s" 1 . Evidently, the increased 
truncation radius results in lower values of the velocity dis- 
persion normalisations. At the 99 per cent confidence level 
we find al < 195 and a s J < 117 km s _1 . 

The results obtained in this section are summarised in 
Tables [2] and El 
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Table 2. Summary of best-fitting models. 



Studied aspect of the model 


Velocity dispersion scaling law 


Truncation radius scaling law 






n 


a* (km s ) 


7 


r* (h 1 kpc) 


x 2 


No lensing 


(7 = 









168.66 


Truncation radius scaling law (first minimum) 


0.33 (fixed) 


114 


-2.8 


170 


164.50 


Truncation radius scaling law (second minimum) 


0.33 (fixed) 


124 


-0.9 


69 


164.66 


Velocity dispersion scaling law 


0.45 


122 


(fixed) 


40 


164.30 


Velocity dispersion scaling law 


0.25 


56 




n = rvir 


168.13 


Star forming/passive velocity dispersion scaling law 


0.33/0.25 (fixed) 


109/196 


(fixed) 


45 


164.59 


Star forming/passive velocity dispersion scaling law 


0.33/0.25 (fixed) 


33/113 




n = r vil 


167.49 



Table 3. Summary of constrains (with confidence level within parenthesis) on velocity dispersion scaling laws after marginalisation. 



Studied aspect of the model Velocity dispersion scaling law 

•q o* (km s —1 ) 

Truncation radius scaling law 0.33 (fixed) 

Velocity dispersion scaling law (7 = 0) 0.45±°' 1% (68.3 p.c.) < 156 (95 p.c.) 

Velocity dispersion scaling law (rt = r v j r ) 0.25^q'29 (68.3 p.c.) < 91 (95 p.c.) 

Star forming velocity dispersion scaling law (7 = 0) 0.33 (fixed) 91^63 (68.3 p.c.) 

Passive velocity dispersion scaling law (7 = 0) 0.25 (fixed) 181^gg (68.3 p.c.) 

Star forming velocity dispersion scaling law (rt = r v j r ) 0.33 (fixed) < 101 (95 p.c.) 

Passive velocity dispersion scaling law (r t = r v ; r ) 0.25 (fixed) 109^1 70 (68.3 p.c.) 
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Figure 5. Magnification-residual diagram for the best-fitting 
halo model with different scaling laws for passive and star forming 
galaxies. Solid squares represent the weighted average residual 
computed for SNe la in magnification bins of width 0.02 mag. An 
average residual error bar is shown in the upper left corner. The 
slope of the solid line is unity. 



5 EXTINCTION BY DUST 

The brightness of SNe la can be affected by dust extinction 
in addition to gravitational lensing. Extinction by dust, in 
contrast to gravitational lensing, leads only to dimming of 
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Figure 6. Best-fitting values (circle) and 68.3, 95, and 99 per cent 
confidence level contours for the velocity dispersion normalization 
for passive (cj) and star forming (crj f ) galaxies after marginali- 
sation over r*. 



SNe la. Extinction could therefore lead to underestimation 
of the effect of gravitational lensin g magnification. 

Recently iMenard et al. I (|2009D used the correlation be- 
tween the brightness of tens of thousands of quasars and the 
position of tens of millions of galaxies in the foreground to 
measure magnification and extinction simultaneously. Ac- 
cording to their measurements the average extinction and 
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magnification have the same dependence on the impact pa- 
rameter (oc £ -0 ' 8 ). 

Dust extinction and gravitational lensing can be sep- 
arated because the former effect is wavelength dependent 
while the latter is achromatic. In order to calibrate the 
SNe la to become standard candles, the _B-band magnitude, 
m_B, is corrected l|Astier et al.ll2006T ). 



[IB = TUB - M + a(s B - 1) - PC, 



(23) 



using stretch, sb, which parameterizes light curve shape, 
and colour, C. The parameters M, a, and /3 are deter- 
mined from the data. The colour correction could poten- 
tia lly correct for some of the du s t ext inction. According 
to IMenard. Kilbinger fc Scranton (2009), who have inves- 
tigated the consequences of dust extinction for SN la cos- 
mology, this colour correction is not sufficient. 

The values of the truncation radius preferred by the 
data are smaller than the virial radius (see section |4}. This 
rather surprising result might be explained by dust extinc- 
tion not accounted for by the colour correction in equa- 
tion (|23[l . Another indication of the small values of rt pre- 
ferred by the data, is the small values of a» and a* f ob- 
tained when assuming r t = r v i r (see the end of section |4~31 . 
We have investigated if extinction by dust could explain the 
small values of these halo parameters. 

First we use the high As model in 
IMenard, Kilbinger fc Scranton! |2009), which only de- 
pends on zsn (see their fig. 1), to estimate the extinction 
in the B-band, Ab, for each SNe la. The values of As are 
then used to correct raj and C [this variable is essentially 
E(B — V)\ assuming a standard dust law. For the high 
A B model Rb = Rv + 1 = 4.9. Thereafter we fit f2 M , M, 
q, and /3 to the extinction corrected data and compute 
new extinction corrected residuals. These parameters 
are hence allowed to compensate for the effect of dust 
extinction. A flat Universe is assumed, i.e. Qa = 1 — Qm- 
In table U the differences in these parameters with re- 
spect to the uncorrected case are listed. Correcting the 
SNe la for redshift dependent extinction has clearly a 
very small effe ct on the parameters, and con s equen tly on 
the residuals. IMenard. Kilbinger fc Scranton! (J2009) also 
find small differences in M, a, and /3 (see their table 1). 
For SIm, however, we find only half of their difference 
(AilM = 0.017). The statistical significance of the difference 
is, however, very similar: they find 0.55a, whereas we fi nd 
0.50a. Note that the parameter /3 ~ 3 l| Tripp! 119981 ) is 
different from Rb = 4 . 1 measured fo r Milky Way dust 
Savage fc Mathid Il979l; ISeatonl Il979l ; iRieke fc Lebofskvl 



ll985l ; ICardelli. Clayton fc Mathislll989h 

Finally, we redo the fit of the halo parameters af and a* 
using the extinction corrected residuals. We find the differ- 
ence in the best-fitting values of the halo parameters due 
to the extinction correction to be negligible. This is not 
very surprising because the difference between the extinc- 
tion corrected and uncorrected residuals is on average only 
0.000 ± 0.005 mag. 

Since we have information about the foreground galax- 
ies, the effect of a more detailed extinction correction, where 
the contribution from each individual galaxy is added to the 
total exti nction, can be studie d. According to the measure- 
ments of IMenard et al] (|2009l ) the mean extinction profile 



around the galaxies in their sample, which span the bright- 
ness range 17 < m; < 21, follow 

-0.84 



(A v )(£) = 4.14 x 10" 



WOh- 1 kpc 



(24) 



The extinction in the B-band can be obtained from the V- 
band extinction via Ab = (Rb / Rv)Av ■ We assume the 
value Rv = 3.1 corresponding to Milky Way dust. For our 
sample of SNe la the more detailed model predicts larger 
extinctions ({Ab) = 0.05 mag) than the redshift dependent 
model ((Ab) = 0.03 mag). Also for this model the differences 
in JIm, M, a, and /3 are small compared to the uncorrected 
case (see table |4|. For the residuals the mean difference is 
0.000 ±0.009 mag. The value of /3 is, however, slightly larger 
than for the redshift dependent extinction model. The halo 
parameters a£ and af are slightly increased by 10 km s _1 
each, when the extinction corrected residuals are used in the 
fitting procedure. 

Consequently, neither of the two dust extinction mod- 
els have a significant impact on the results. We conclude 
therfore that dust extinction is not a plausible explain for 
the small values of the velocity dispersion normalisations we 
find for r t = r v i r . 



6 GRAVITATIONAL LENSING BRIGHTNESS 
SCATTER 

The results of the previous sections allow us to investigate 
the implications for supernova cosmology. We restrict our- 
selves to the model investigated in section l4~3l where passive 
and star forming galaxies are described by different velocity 
dispersion scaling laws. This model is more realistic than 
the ones where all galaxies are treated as if they were of the 
same type. 

As already stated in the introduction, gravitational 
magnification leads to additional SN la brightness scatter. 
This scatter is expected to increase with redshift. The thick 
solid curve in Fig. [7] indicates the dispersion in Ami cns as a 
function of redshift computed for the best-fitting model af- 
ter marginalisation over r* using randomly selected lines of 
sight. We use the best-fitting truncation radius, r» = 45/i _1 
kpc. The light and dark shaded regions indicate the scatter 
at the 68.3 and 95 per cent confidence level translated from 
Fig. El 

Using simulations iHolz fc Linderl l|2005h predicted the 
increase in lensing dispersion to be approximately propor- 
tional to the SN la redshift, <7Am lcns ~ Bzsn, with B = 0.088 
mag. In Fig.[7]we show the predicted dispersion for constants 
of proportionality in the range B — 0.01, 0.02, . . . , 0.15 mag. 
From the figure we conclude that it is reasonable to de- 
scribe the redshift dependence of the lensing scatter as pro- 
portional to the SN la redshift. Ot her authors hav e foun d 
lower values for the di spersion than lHolz fc Linderl |2005). 
IBergstrom et all (|2000h find ~ 0.04 mag for zsn = 1 

using simulations and assuming smooth halo profiles. 



From Fig. [3 we find B ~ 0.055] 



mag. At the 95 per 



cent confidence level we find B < 0.1 2 mag. The results are 
thu s consistent with the sim ulations of lHolz fc Linderl (I2005T ) 
and IBergstrom et al.l (|2000h at the 68.3 per cent confidence 
level. 

The value of the gravitational lensing dispersion is very 
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Table 4. Differences in cosmological parameters due to dust extinction correction compared to the uncorrected case. 
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Extinction correction 


Aa 


A/3 


AM 


AQ M 


Redshift dependent 


0.002 (0.03ct) 


0.003 (0.03cr) 


0.001 (0.08a) 


0.009 (0.50a) 


From individual galaxies 


0.003 (0.04tr) 


0.026 (0.26ct) 


-0.001 (0.08a) 


0.007 (0.39a) 



sensitive to the high magnification tail. For small samples 
of SNe la, which fail to sample the high magnification tail, 
we might therefore find a smaller dispersion. This appears 
to be the case for our sample of 175 SNe la. The circles in 
Fig. [7] shows the lensing scatter of these SNe la computed 
for four redshift bins with ~ 44 objects in each bin. The 
circles are consistent with B ~ 0.035 mag, which is smaller 
than the value predicted by the best-fitting model. Reject- 
ing simulated SNe la with Ami cns < —0.25 mag leads to the 
dispersion shown by the thick dashed curve, which agrees 
with the scatter for the SN la sample. Only 0.3 (0.7) per 
cent of the SNe la would be brighter than —0.25 mag for 
SNe with redshift 0.7 (1.0). For the whole sample of 243 
SNe la the best-fitting model predicts the number of super- 
novae brighter than —0.25 mag to be 0.8. For the subsample 
of SNe la used in the lensing analysis the corresponding 
number is 0.6. Since the number of SNe la belonging to the 
tail is expected to be very small, our sample of SNe la is 
probably not sampling the high redshift tail brighter than 
Amiens — —0.25 mag. 



7 DISCUSSION AND SUMMARY 

We have used presumably lensed SNe la from 3-year of 
the SNLS to investigate properties of dark matter haloes of 
galaxies in the deep CFHTLS fields. The dark matter haloes 
were modeled as truncated singular isothermal spheres with 
velocity dispersion and truncation radius given by luminos- 
ity dependent scaling laws. 



Another way to probe da rk matter haloes 



is 


via galaxy-galaxy lensing (ITvson et al.l 


19841; 


Brainerd. Blandford & Smail 19961; Hudson et al.l 


1991; 


Fischer et al.l l200d: iMcKav et all l200ll; iGuzik & ! 


Seljak 


2002 


; Hoekstra et al. 20031. 20041; iKleinheinrich et al. 


2006 


; IMandelbaum et al.l 120061: iParker et al.1 120071; 


Mandelbaum, Uros & Hiratal 120081; iMandelbaum et al. 


2009T). which relies UDon measurement of shear via the 



ellipticity of lensed galaxies rather than the convergence. 
In the galaxy-galaxy lensing literature different truncation 
radius scaling laws o f the form expr e ssed b y equation (|15[) 
have been explored. iHoekstra et al.l (|2003T ) considered, for 
example, the size of galaxy haloes to be either the same for 
all galaxies (7 = 0) or scale as a 2 (7 = 2). Since the virial 
radius of a SIS is proportional to its velocity dispersion, 
another plausible scaling law would be 7 = 1. 

Fixing the exponent of the velocity dispersion scaling 
law, equation ()13p . to 77 = 1/3 we tried to constrain the 
truncation radius scaling law, but the data did not allow us 
to constrain r* or 7. 

We have also explored the velocity dispersion scaling 
law b y setting 7 = 0. Figure [3] shows the result of this ex- 
ercise. IHoekstra etHl (|2004T ) used galaxy-galaxy lensing to 
investigate the properties of dark matter. They used a model 



ID 




0.2 0.4 0.6 0.8 1 

2 SN 

Figure 7. Gravitational lensing scatter as a function of red- 
shift. The thick solid curve shows the dispersion for the best- 
fitting model found in section 14.31 Light and dark shaded re- 
gions correspond to 68.3 and 95 per cent confidence level, re- 
spectively. Predictions for a simple model of the lensing scatter, 
aAm lenB = Bzqn, for B = 0.01,0.02, . . . ,0.15 mag are indicated 
by the thin lines. The circles show the scatter in the magnification 
computed in four redshift bins for the sample of SNe la used in 
the analysis. Horizontal error bars indicate the width of the red- 
shift bins. The thick dashed curve shows the dispersion for the 
best-fitting model when SNe la with Attisn < —0.25 mag have 
been discarded. 



slightly differen t from ours, namely the truncated iso thermal 
sphere model ferainerd. Blandford fc Smaill Il996h . which 
has a parameter s describing the truncation scale of the halo. 
For a value of r\ fixed to 0.3 they found a, — 136 ± 5 ± 3 
km s" 1 (statisti cal and systematic uncertainties) and s = 
185128 k P c - IParker et all (|2007T ) found a similar value 
of the velocity dispersion, a* = 137 ± 11 km s , us- 
ing CFHTLS data. I n ano ther galaxy-galaxy lensing study 
IKleinheinrich et all |2006l ) found <r» = 132^4 km s" 1 and 
r/ = 0.37±0.15 for SIS haloes truncated at 350ft" 1 kpc. Our 
results for the velocity dispersion scaling law are in good 
agreement with the findings of these studies. The data pre- 
fer, however, a much smaller truncation radius. 

Moreover, we tried to disentangle the scaling laws of 
passive and star forming galaxies. Again we restricted our- 
selves to models with 7 = 0. Passive and star forming galax- 
ies were assumed to obey F-J and T-F relations, respec- 
tively. The best-fitting values of the velocity dispersion and 
confidence level contours obtained after marginalisation over 
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r» are shown in Fig. [6] The point <r» = af = in this plot, 
which corresponds to no lensing, is excluded at the 91.6 per 
cent confidence level. 

The results agree well with empirical F-J and T-F rela- 
tions. iMitchell et al] (|2005h measured the F-J relation in the 
r-band using Sloan Digital Sky Survey data. Their measure- 
ment corresponds to a velocity dispersion normalisation of 
190 km s _1 fo r our choice of L* , i .e. Mb = —20.3 mag in the 
Vega system. Bohm et al.l (|2004h de rived a redshift depen - 
dent T-F relation using data from IPierce fc Tullvl l|l992h . 
For a galaxy at z = their findings corresponds to a nor- 
malisation of 122 km s _1 . 

Assuming the truncation radius to be given by the 
virial radius we found al < 195 and af < 117 km s _1 , 
a result which is marginally consistent with haloes trun- 
cated at the virial radius with velocity dispersion given 
by F-J and T-F relations. Corr ecting the SNe la for dust 
extinction using the results o f Menard et al.l |2009l ) and 
iMenard. Kilbinger fc Scrantonl (2009) have no significant 
impact on al and <r» f . 

Furthermore, we have used the best-fitting model with 
different scaling laws for passive and star forming galax- 
ies to investigate the gravitational lensing scatter. For a 
model where the scatter is proportional to source redshift, 



Bzsn, we find B ~ 0.055 



+0.039 



mag (at 68.3 per 



0"Am le 

cent confidence level and B < 0.12 mag at the 95 per cent 
confide nce level). The scatter is consiste nt wit h the predic- 
tion of iHolz fe Linderl (B = 0.088 mag. l2005h at the 68.3 
per cent confidence level. 

If the best-fitting model gives an accurate description 
of the gravitational lensing of SNe la, the contribution 
from gravitational lensing scatter to the intrinsic dispersion 
(~ 0.09 mag) of the SNLS SNe la is very small. For SN la 
data sets with zsn > 1.6, gravitational lensing could, ac- 
cording to the best-fitting model, contribute significantly to 
the brightness scatter. The farthest known SN la (1997ff) 
with zsn — 1.8 is indeed believed to be magnified by a 
few t e nths of a magnitude llLewis fe Iba ta 2002; Ricss et al 



2001 



2002 



Mortsell. Gunnarsson fe GoobarT 200 ll ; Bemtez et al 
I Jonsson et al.ll2006l ). 
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